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ABSTRACT 

Laboratory measurements show that the spectra of the capillary 
waves grow with wind speed over six orders of magnitude. Measure- 
ments in a wind-water tunnel show a simple one to one relationship be- 
tween cryy and wind speed. The laboratory data for both waves and 
backscatter have negligible scatter because the wind for each measure- 
ment was constant and the radar return and wave spectra were averaged 
over a long enough time. Other conditions did not change. 

Regression analyses of upwind and crosswind values of 
10 logiQ (a°(9)/rr*°(10 ° )) from aircraft data show linear relationships 
with trie logarithm of the wind speed at 15°, 25° and 35°. The data have 
substantial scatter in part because both the average wind speed and 
average sea return values had to be estimated from relatively short 
samples. The scatter can be partially understood and predicted from 
a combination of turbulence theory, radar theory, and the small sampl- 
ing theory of statistical inference. The scatter is caused in part by 
two effects; the turbulence of the wind over the water and the ’’fading” 
characteristics of the radar return. 

When these results are applied to a prediction of the sea return 
values to be obtained by S193 on Skylab, it can be shown that the size 
of the illuminated patch effectively averages out the horizontal scales 
of gustiness so that the measurement will correspond to the synoptic 
scale wind. The effect of ’’fading 11 will introduce a scatter of less than 
six percent for angles greater than 15°. 

There will still be some difficult problems to be solved for S193. 
These problems are described, and methods for their solution are re- 
viewed. 
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INTRODUCTION 


Both laboratory and aircraft measurements over the ocean increas- 
ingly support a strong correlation between radar sea return and wind 
velocity. It is dangerous to oversimplify the results obtained in the 
laboratory, and apply them directly to the open ocean without first ac- 
counting for differences between laboratory and open ocean conditions. 
However, once these differences are accounted for, the combined in- 
formation obtained from the laboratory and from aircraft measurements 
makes it possible to understand some of the reasons for the scatter in 
the aircraft data and to calculate what these effects will be in the mea- 
surements to be made by SI 93 on Skylab. Such an analysis also reveals 
that only part of the scatter in the aircraft data is caused by sampling 
variability in the data. Reasons for this additional scatter are postu- 
lated. Many of these other sources of scatter will not be present in 
SI 93, and still others can be removed with the larger data base to be 
made available when Skylab is flown. 


LABORATORY DATA 


Laboratory data on the form of the capillary wave spectrum have 
been combined to show that the capillary spectrum in a wind -water tun- 
nel varies over six orders of magnitude as shown by Pierson, Jackson, 
et al. (1971). In Figure 1, the capillary waves are very low below a 
friction velocity, u$, of 10 or 12 cm/sec. Then the waves grow by 
more than four orders of magnitude as the friction velocity changes by 
just a few cm/sec. The spectrum still grows with increasing wind 
speed all the way to friction velocities of 150 cm/sec, but more slowly. 
Photographic data also support the conclusion that the spectrum grows 
with increasing wind speed in the capillary frequency range. 

Wright has measured &yy at 60° off the vertical in a wind-water 
tunnel for a wide range of wind speeds. These values are the exact 
unnormalized values for the scattering cross section. Some of the 
preliminary results are shown in Figure 2. The backscatter from the 
capillary waves clearly increases with increasing wind speed. The 
change from a very steep slope to a less steep slope in this plot corres- 
ponds to the value D equals one in Figure 1. Wind speeds in the free 
air section of the tunnel are given in feet/min. An auxiliary scale gives 
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FIG. 2 RADAR BA C KSC AT TE R AS A F UN CT ION OF 
WIND SPEED IN A WIND-WATER TUNNEL. DATA 
COURTESY OF DR. JOHN WRIGHT. 
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the approximate equivalent wind at 1 0 meters over the ocean. 

In experiments in wind -water tunnels, it is possible to keep condi- 
tions relatively simple compared to the open ocean. Measurements 
can be averaged over a sufficiently long time so as to ensure the elimi- 
nation of sampling variability. Near the surface of the open ocean, wind 
speeds of the values shown would produce much larger waves, with a 
spectrum extending to much lower frequencies. Moreover, the winds 
over the ocean are not as steady as the winds in a wind-water tunnel. 


ATMOSPHERIC TURBULENCE 


The atmosphere is turbulent; its many scales of motion can be as- 
signed to three categories; the microscale, the mesoscale and the 
macroscale, or the planetary scale. A wind vane and an anemometer 
respond to the frequencies in the microscale with some damping. The 
mesoscale is often unimportant over the oceans, but it is important for 
the study of, for example, mountain- valley winds, thunderstorms, and 
heat islands over cities. 


Lumley and Panofsky (1964) cite an important paper by Van der 
Hoven (1957) that first attempted to study the full range of motions for 
the horizontal wind. A figure from Lumley and Panofsky (1964) which 
is redrawn from Van der Hoven (1957) is repeated here as Figure 3. 


This figure is an oversimplification of the problem; yet the over- 
simplification provides an important insight into the nature of turbu- 
lence. The vertical scale is nS(n) where n is the frequency and S(n) 
is the variance spectrum of the wind with dimensions of (velocity)^ 
times time so that the vertical axis has the dimensions of m^ sec"^. 
The horizontal scale is a logarithmic scale in n where n is in cycles 
per hour. 


The graph has the property that 

n„ n„ 

2 2 

^ nS(n) d log n = ^ n S(n) 
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so that this form preserves area in contrast to one where the logarithm 
of S(n) is plotted against the logarithm of n. 

Figure 3 has three relative maxima; one near 50 cycles per hour, 
one near one cycle per twelve hours and one near one cycle per four 
days. The part on the left is typical of what might be computed for S(n) 
if say, the winds were recorded by an anemometer near the ground once 
an hour as one hour averages for a year or so and then the variance 
spectrum was computed. The part on the left would vary considerably 
from one place in the world to another. 

The part on the right represents the spectrum of the gusts and lulls 
of the horizontal wind. It really depends to a very large extent on the 
speed and direction of the average planetary wind during the time the re- 
cord was obtained from which the spectrum on the right side was com- 
puted. If the synoptic scale wind is high, the spectrum on the right will 
be high. If the wind is low, the spectrum will be low. 

The usual anemometer installation is hardly the optimum instrumen- 
tation for the study of turbulence, but studies of turbulence with faster 
response instrumentation by Phelps and Pond (1971) and DeLeonibus (1971) 
show that simple anemometers yield useful information on the important 
scales of turbulence in the wind. Much of what was known through 1964 
is summarized by Lumley and Panofsky (1964). 

The important point about Figure 3 is not the three peaks in the spec- 
trum, which fluctuate in amplitude and position, but the low region in the 
spectrum extending from about .05 to 10 cycles per hour. The area under 
the curve to the left of this band is five to six times greater than the area 
under the curve in the band. The area to the right is four to five times 
greater . 

The low region in this spectrum justifies the fundamental hypothesis 
of numerical weather prediction; namely that the turbulent high frequency 
part of atmospheric motions can be filtered out of the synoptic scales, and 
parameterized as "eddy Austausch" coefficients so that numerical weather 
predictions can be carried out on a grid with spacings large compared to 
the scales of the eddies in the spectra on the right hand side. The numeri- 
cal weather prediction schemes use functions of the synoptic scale winds 
that depend on a parameterization of the turbulence. The synoptic scale 
winds are nearly horizontal with the vertical components of the motion 
being 10“^ times the horizontal components. 
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If V(x,y, z, t) is the horizontal wind in the planetary boundary layer 
over the ocean, then 


— t+T/2 

V(x,y, x, t) = Tp, ^ V(x,y,z, t ) dt (2) 

1 J o o 

t-T/2 

and 

V (x, y, z, t) = V (x, y, z, t) + V* (x, y, z, t) (3) 

where 

V 1 (x, y, z, t) = V(x, y, z, t) - V (x, y, z, t) (4) 

The turbulence, that is the fluctuations in the wind to the right of, 
say, one cycle per half hour, is now in the V . The synoptic scale wind, 
in the absence of any important mesoscale perturbation, such as thunder- 
storms, is now the V. 

Equally well, one could write equation 5 instead of equation 2 and 
substitute a space average for a time average where A is some area 
over the ocean surface, 

V (x, y, z, t) =J- ^ V( X o’ y Q ’ z* t) dx Q dy Q (5) 

A 


These concepts have been applied to a study of the best way to pro- 
cess anemometer data to be obtained by the National Data Buoy System 
of NOAA by Baer and Withee (1972a, 1972b). Numerous spectra similar 
to the one shown in Figure 3 for Argus Island, Keflavik, Iceland, and a 
site in the Pacific Ocean are shown that have the same properties as 
Figure 3. A report by Millard (1968) is cited that also show a similar 
spectrum. 

The problem addressed in this study is to provide the best estimate 
of the horizontal vector wind for numerical forecasting purposes and for 
synoptic meteorological purposes. It is shown that the best that can be 
done is to give the wind, for a homogeneous area of flow, to within a 
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standard deviation of about one knot. 

Equation 2 represents the simplest way to estimate a synoptic scale 
wind. It corresponds to the so-called ir box-car" or rectangular weight 
function in which the limits of integration are set from -« to +o> and a 
function G(t - t Q ) equal to 1 /T for t - T /2 < t Q < t + T /2 is inserted 
after the integral sign. This weight function is not too good in its effect 
on the spectrum of the wind because it does not act like a sharp low pass 
filter. Frequencies both below 2tt/T and above 2tt/T are passed in 
varying amounts by the filter. A properly normalized weight function of 
tne form cos Tr(t- t 0 )/T can be shown to improve the filter characteris- 
tics, for example. 

It is to be expected that improved ways to process anemometer data 
will result from further work in this problem area. For example, Baer 
and Withee (1972b), have studied how to locate and report the sudden 
changes, as a function of time, in wind direction and speed that occur 
during frontal passages. 

For the purposes of this study of the aircraft program and of the re- 
sults to be expected from Si 93, it should be noted first of all, that a 
scattering of data buoys over the ocean comparable to the scattering of 
ships in the results of Druyan(197l) could provide calibration points for 
a radar -radiometer system on a spacecraft in the form of properly 
averaged vector winds good to about 1 knot standard deviation in speed 
and a few degrees in direction. 

It should also be noted that there is room for improvement in the 
area average indicated by equation (5). The average need not be re- 
stricted in concept to an average of weight one inside a certain area and 
zero outside the area. A properly chosen weight function can filter the 
wave number spectrum of the horizontal wind for numerical prediction 
purposes very efficiently. 

Taylor's Hypothesis . - Taylor's hypothesis provides a means to 
relate measurements of the wind at a point as a function of time to mea- 
surements of the wind along a line (usually in the upwind downwind direc- 
tion) as a function of distance. One way of describing the hypothesis is 
to say that it assumes that the pattern of turbulence is"frozen n and is 
advected past a point by the mean wind. This is equivalent to the as- 
sumption that the length of an eddy in space is given by multiplying the 
period of an eddy in time by the mean wind speed. 
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Lumley and Panofsky (1964) state that "Lappe, Davidson and Notess 
(1959) show by a comparison between tower measurements and airplane 
measurements at 90 and 120 m high, respectively, (the higher speed air- 
plane measurements are almost sure to be satisfactory) that the hypo- 
thesis appears to be valid for horizontal wind fluctuation but that the 
vertical fluctuation measured by airplanes contains more low-frequency 
energy than the tower measurements. The authors ascribe this discrep- 
ancy to the effect of wind shear described by Lin. On the other hand, 
however, the terrain is not homogeneous, and flow patterns fixed rela- 
tive to the ground contributing to the airplane, but not to tower measure- 
ments, may be responsible for the discrepancy". 


"In summary, the very limited information seems to indicate that 
Taylor’s hypothesis is generally applicable except for the low wave 
number components of vertical velocity fluctuations. Therefore all time 
spectra and correlation functions observed at a fixed point or from air- 
planes will, in what follows, be interpreted as space spectra and cor- 
relation functions along the direction of mean motion". 

More recently, multiple arrays of wind sensors seem to suggest that 
eddies are dispersive. Typically, eddies can travel about six times 
their own length before they die out as shown by McDonald and Herrin 
(1971). The use of the Taylor’s hypothesis in what follows is conserva- 
tive. The additional decorrelation that is observed would increase the 
estimates of the number of degrees of freedom. 

Argus Island Anemometer Data. - Figure 4 shows a small segment 
of an anemometer record obtained at Argus Island during Mission 119 
(1970). This particular segment is a portion of the record obtained dur- 
ing the time near Flight 7, Run 3, Line 1. It shows the winds fluctuating 
from 12 knots to 21 knots in 60 seconds. A total of 2000 values of this 
Argus Island anemometer record were read every 4 seconds 1000 values 
before the Flight 7, Run 1, Line 3 measurements and for 1000 values 
after . 

The covariance function for the anemometer record is shown in 
Figure 5. The variance spectrum is shown in Figure 6 as plotted on a 
relative scale. The covariance function shows a quasi-exponential non- 
oscillatory decay out to 144 seconds as plotted. At 400 seconds, the 
covariance is still equal to 0. 28 or about 7$ of the value at the origin. 
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Each spectral estimate has about 40 degrees of freedom according 
to Blackman and Tukey (1958). The spectrum shows that the major por- 
tion of the variability is concentrated in the low frequencies, corres- 
ponding to periods of 50 seconds and higher. 


A well-known problem in statistics, is that of estimating the mean 
and variance of a normal population given a sample of N independent 
observations from a given population and of determining the variance 
and standard deviation of the estimate of the mean. For a sample of 
N independent observations the variance of the sample mean is given 
by the estimate of the variance (SD V )^ of the population divided by N 
as in equation 6. 

. (SD r 


cr = 
V 


N 


( 6 ) 


The 2000 points in the anemometer record are not independent and 
therefore equation 6 cannot be used with N equal to 2000. There are 
several ways to account for the correlation between points in such a 
time series. An equation given in Neumann and Pierson (1966) as given 
originally in Blackman and Tukey (1958) accounts for the decorrelations 
at different lags and provides for an effective number of degrees of 
freedom N v as in 


n" 


N 

m 


[SS*1 2 

vi, 2 

2(S V ) 


(7) 


where m is the number of smoothed spectral estimates available and 
the S* are the hanned spectral estimates. This equation holds best 
for spectra that do not depart too much from a white noise spectrum. 

For this particular anemometer record, the effective number of 
degrees of freedom may be only about 178. The 2000 point autocorre- 
lated record is perhaps the equivalent of a totally uncorrelated sample 
of 178 values from a normal population. The standard deviation of the 
sample mean is thus quite a bit larger than it would be for the same 
variance if the sample had consisted of 2000 independent points. 

If the 178 independent points were thought to be equally spaced over 
the 2000 point time history, they would be located at about every 1 1th 
point, or every 44 seconds. This statement needs to be interpreted 
properly. Clearly from Figure 5 the covariance is still fairly high at 
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44 seconds. The real point is that this complicated, non-independent 
time series of points contains enough information on the variance of the 
process and its frequency content to make it possible to replace it by an 

•A. 

equivalent NT independent points for the purpose of determining con- 
fidence intervals on V. The 44 second apart effectively independent 
points in the anemometer record, if advected at a speed of 8 meters per 
second, implies that points 350 meters apart along an upwind line in 
space would be effectively independent in this same sense. 

General Properties of the Turbulent Wind. - It is just as difficult 
to get a reliable estimate of the mean wind when a radar measurement 
is made over the ocean as it is to get a reliable estimate of the radar 
scattering cross section. In the plots and tables to follow, the scatter 
is caused just as much by inadequacies in the wind speed measurements 
as it is caused by the possible errors in the radar measurements. 

In working with the turbulent wind, meteorologists study both the 
spectra and cross spectra of the V' components as described above 
and the vertical velocity fluctuations, w 1 . They also study the varia- 
tion of the mean wind as a function of elevation above the ocean surface, 
V(z). 


The variation of mean wind with height above the ocean surface de- 
pends on the stability of the layer of air over the water, which in turn 
is controlled by the air -sea temperature difference. If the air -sea tem- 
perature difference is small, the lapse rate in the atmosphere is adia- 
batic and the wind varies with height according a logarithmic wind pro- 
file. If the air is colder than the water the atmosphere is unstable, and, 
for the same wind above the surface there will be a greater stress on 
the surface. Conversely, with the air warmer than the water the atmos- 
phere is stable and for a given wind above the surface, it will be less. 

All of these factors were taken into account by Cardone (1969) in 
developing a computer based procedure for studying the winds in the 
planetary boundary layer. They are also taken into account in the inter- 
pretation of the winds for the various missions. The winds reported 
and used in these studies are the ’'equivalent 11 winds for a neutrally 
stratified atmosphere referred to 19.5 meters above the ocean surface. 
They would produce the same stress at the surface as the measured 
wind at some other elevation with some known air sea temperature dif- 
ference according to the available theories on the subject. In general, 
the effect of air -sea temperature difference is relatively unimportant 
for high winds. 
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The fluctuations in wind speed recorded by a cup anemometer are 
more or less the equivalent of the fluctuations of the vector component 
of the wind in the direction of the mean wind. The spectrum shown in 
Figure 6 is typical of such spectra in that a record 8000 seconds long does 
not resolve the maximum away from the origin. 

In general, according to Lumley and Panofsky (1964), the variance 
of this turbulent spectrum grows as the square of the mean wind as in 
equation 8 

00 

VAR = (SD v ) 2 = ^ S(n) dn a: K V 2 ( 8 ) 

o 

when the wind is measured at a fixed height. 

The value of SD V is relatively constant with height so that, since 
the mean wind increases with height, the "relative gustiness" SD v /V de- 
creases with height. For a fixed height SD v /V is_also dependent on the 
air-sea temperature difference. In general, SD v /V can be considered 
to range from about 0. 10 to 0. 30 with high values reserved for unstable 
air and light winds. 

For the anemometer record being studied the mean wind was 17.0 
knots before correction for height above the surface and stability effects 
and the standard deviation was 1.95 knots so that, at anemometer height, 

SD 

77 ^ =0.114 ( 9 ) 

V 


THE AIRCRAFT PROGRAM 


Radar Sea Return Measurements and Wind Speed. - The scattero- 
meter used to measure sea return was a continuous wave radar with a 
vertically polarized fan beam that differentiated the angle 9 off the 
vertical by means of Doppler shifts in the backscattered radiation. 
Bradley (1971) has carried out an extensive analysis of this system. 
This study relies heavily upon his analysis, differing however, in some 
important points. 
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A typical upwind data run would be obtained by flying for about five 
minutes in a direction chosen visually to be perpendicular to the direc- 
tion of the oncoming waves at the sea surface. The backscatter ed signal 
from ahead is upshifted in frequency, while that from behind is down- 
shifted. Appropriate signal processing techniques yield measurements 
of 0*0 in the vertical transmit vertical receive mode for various angles 
off the vertical. 

The largest slice of recorded data that could be handled by the avail- 
able computer at MSC was 0. 3275 seconds long, and every fourth one 
of these data slices (spaced consequently 1.31 seconds apart) was chosen 
for analysis. Each 0.3275 second data slice yielded a measurement of 
c r° at all desired angles. A four minute data run would yield about 180 
separate measurements of cr° for each angle. For such a run, a single 
average value of o*° obtained from the individual values computed from 
each 0. 3275 second sample was obtained for comparison with the average 
measured, or estimated, wind speed. Data processed for the different 
runs varied in duration, so the number 180 is appropriate only for a 4 
minute run. 

_Bradley (1971) has given the result of a regression analysis of cr° 
on V for the data obtained during Missions 119 and 156. 

The values of A and B in equation (10) were determined for upwind 
and cross wind conditions for 35°, 25°, and 15° as a least square re- 
gression of log V on a^. The values of A and B are given in Table 1. 

10 lo g in = A+B log v 00) 

1U tr° (10°) 1U 


TABLE I.- RESULTS OF CURVE FITTING BY 

BRADLEY (1971) FOR MISSIONS 119 AND 156 


Wind Direction 

e 

A 

B 

RMS Error 

upwind 

15 • 

- 7.51 

3. 30 

0. 717 db 


25 ° 

-25.56 

11.21 

0.892 


35 ° 

-34. 42 

14. 93 

0.882 

crosswind 

15° 

- 6.68 

2. 1 3 

0.748 


25 ° 

-22.84 

7. 51 

1.196 


35 ° 

-35.51 

12.62 

1.355 
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The use of equation 10 to study the individual values of <r° provides 
a great amount of information on the nature of the return. Denote 
o-°(9)/o-°(l0 0 ) by eg (9) in equation 10 so that 


10 log l0 (cr N (0)) = A +BlDg 10 V 


(ID 


A change of in ojsj must produce a corresponding change of AV 

in Y as in equation 12 . 


, o 


o . 


10 lo g ]0 (<r N + A<t n ) = A + B log 1() (V + AV) 


— ' AV 

= A +B log i() V + B log ]() (l + 37 j (12) 


If equation 11 is subtracted from equation 12, the result is equation 13. 


A«r, 


°N 


B/10 


(13) 


Changes in (one of the individual values observed as a time 

series during a particular flight line) are caused by two effects: ( 1 ) 
the radar backscatter is a random variable because of the noiselike 
character of the fading and ( 2 ) the backscatter ed signal strength changes 
because of actual changes in the roughness of the surface. 

The variability of the measurements of <r° due to fading depends on 
the total number of independent samples in the return signal. Bradley 
(1971) has given the formula for this source of variability and from it, 
it is possible to compute SD^ , the standard deviation of * 

From equation 13, 1 +AV/V and 1 - AV/V can then be computed 
as in Table II. 
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TABLE II.- VARIATION OF APPARENT WIND SPEED DUE TO 
FADING FOR A 0. 3275 SECOND MEASUREMENT 
FOR AN UNDERLYING SURFACE WITH CONSTANT 
ROUGHNESS. 


CD 

1+ ti ! 

°n 

1 + AV 
V 

1 - AV 
V 

15 e 

1.12 

1 . 41 

0. 71 

25 * 

1.13 

1.12 

0.89 

35 0 

1 . 15 

1.10 

i 

0.91 


The equivalent values in db for 1 + Ao"° /cr in Table II are 0. 49, 0.53 
and 0.62 db. N N 

Fading for this particular instrument over a uniformly rough surface 
amounting to 12^ thus causes 41$ apparent wind fluctuations at 0 equal 
to 15°. At 35°, in contrast, fading produces only 10$ apparent fluctua- 
tions in the wind. For example, if the true wind had been 30 knots, a 
value of cr 0 at 15 ° could have been obtained that would indicate a wind 
as high as 42 knots or as low as 21 knots for about one-third of the values. 
At 35 °, the wind could be computed to be as high as 33 knots and as low 
as 2 7 knots . 

The scattering cross section will also vary if the roughness of the 
underlying surface varies. This roughness is dominantly controlled 
by the capillary wave structure, and as shown in Figure 1, the rough- 
ness increases with increasing wind speed. The capillary waves re- 
spond very quickly to variations in the wind speed so that if the wind 
is, say, 20$ higher than the mean over a particular illuminated cell, 

( r° will increase accordingly. The amount of this change is shown in 
Table III. Possible values for the ratio of the standard deviation of 
the wind speed to the mean wind speed are listed across the top of the 
table. The body of the table contains the resulting change of 
1 + A /( r^ and the equivalent value in db. 
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TABLE III. - VARIATION OF 1 + Ao-°/cr 0 IN TERMS OF 

FLUCTUATIONS OF THE WIND ABOUT THE _ 

MEAN WIND SPEED AS EXPRESSED BY SD /V 

_ v 

SD /V 
v 



0. 15 


0. 20 


0. 30 



, , —o.— o 

1 + A<r /<r 

db 

1 +A<r /cr 

db 

1 +A<r /<r 


15° 

1 . 05 

0. 20 

1 . 06 

0.26 

1.09 

0. 38 

25 ° 

1 . 18 

0. 73 

1 . 23 

0. 89 

1 . 34 

1.28 

35 ° 

1 . 23 

0.91 

1 . 31 

1 . 18 

1 . 48 

1.70 


For example, if SD v /V is 0.20 and the mean wind is say 30 knots, 
there will be patches over the ocean where the wind gusts will be 36 
knots, and at 35°, <r° would increase by 1 . ] 8 db if the average value 
of <r° for a 36 knot wind were recorded. 

Of course, these two effects occur together. A high gust will pro- 
duce sea return that scatters about the expected value for that high wind 
according to the values in Table II. For example, a 30 knot average 
wind could have a local gust of 36 knots where the expected sea return 
would be 1. 18 db higher because of the gust. However, sampling vari- 
ability due to fading would produce fluctuations of ±0.62 db about this 
value . 

It is important to note that the fading effect is large for small 0 
and small for large 9 whereas the variation in wind speed produces a 
large effect for large 9 and a small effect for small 9. 

For future study, the bar over (r° (0) and over V can be removed 
and equation 10 can be written as equation 14. 


10 log 


10 


q ~°(9 ) 

?°(io-) 


= A + B log V 
e l 0 


= 10 log a + B log,. 

e l0 s 10 


V 


(14) 
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Now, for example <r°(0) represents the individual values determined by 
0. 3275 second averages obtained every 1. 31 seconds during an upwind 
flight at, say, 35° and V represents the wind speed that would be com- 
puted by solving for V in equation 14 as in equation 15. 


V (cr° (0 ) ) = 10 



r°(0) 

10 a f (10°) 



o-°(9) 1 

F(IO-) J 


10 

B 


(15) 


Segments of the graphs that result from treating the successive 
values of c r° as a time series and computing V(cr o (0)) for 35°, 25°, and 
15“ for Flight 7, Run 1, Line 1 (an upwind run) are shown in Figure 7. 

The aircraft ground speed should be added to the mean wind speed 
to obtain the speed of the aircraft relative to the "frozen" field of turbu- 
lence. The measurements of <r° , 1 . 31 seconds apart are effectively 

127 meters apart in the turbulence pattern. Were the pattern advected 
by the mean wind at, say, 8 meters per second, equivalent measure- 
ments by an anemometer would be about 16 seconds apart. A four minute 
long aircraft flight is thus the equivalent of a (16/1.31) X 4— 49 minute 
anemometer record. 

These three functions do not look very much like the anemometer 
record of the wind shown in Figure 4. The one for 35 0 looks the most 
like an anemometer record, and the one for 15“ has fluctuations in it 
for an excess of what could be expected in an anemometer record. 

The analysis of these three space histories (after the use of Taylor s 
hypothesis) is shown in the next six figures. Figures 8 and 9 show the 
covariance function and the spectrum for V(<r°(35*)). Figures 10 and 
11 show these functions for 25°, and Figures 12 and 1 3, for 15°. The 
spectral estimates as computed from formulas given by Blackman and 
Tukey (1958) have only seven degrees of freedom so that the confidence 
intervals are quite broad as shown in the figures. 

The curves for 35 “ show that this record has a standard deviation 
of 2.75 knots and therefore a variance of 7.56 (knots)^. The covariance 
function in Figure 8 drops off to about one-half the value at zero at the 
first lag. It does not really get to zero until just short of 1 kilometer. 
The spectrum has a large contribution to the total variance in the wave 
number range below 1 /2036 m - ^ and substantial contributions on to 
1 /508 m~ * . 




0.325 SEC., SAMPLES EVERY 1.31 









97-26 


Figures 8 and 9 should be compared with Figures 5 and 6 from the 
anemometer record. The variance of the anemometer record was 
3.82 (knots) 2 , and the variance of V(<r°(35 °)) is 7.56 (knots) 2 . The re- 
cord graphed as V (cr (35 ')) is thus a combination of the true fluctuations 
of the wind over the water and fading as described by Tables II and III. 

The peak at the left^represents the anemometer fluctuations, which con- 
tribute 3. 82 (knots) to the variance, and the nearly flat part on the right, 
represents the superimposed white noise effect of the uncorrelated fad- 
ing fluctuations which contribute 3.74 (knots) 2 to the variance. They 
happen to be approximately equal for this particular set of circumstances. 

For 25°, the corresponding figures are Figures 10 and 11. The 
variance is 10.96 (knots) 2 so that 7.12 (knots) 2 are due to radar white 
noise effects (fading) and the rest to actual wind fluctuations. The rela- 
tive peak in the spectrum at low frequencies despite the fact that almost 
70$ of the variability in the series is uncorrelated noise is notable. 

For 15°, the corresponding figures are Figures 12 and 13. The 
variance is 110 (knots) 2 and less than 3$ of the variability would be 
attributable to actual fluctuations in the wind over the water. The covari- 
ance function oscillates in a random way above and below zero, and the 
spectrum, since each point has only seven degrees of freedom is essenti- 
ally a white noise spectrum. 

The most important feature of the preceeding analysis is that the 
radar sea return has been shown to be responsive to one, two and three 
knot fluctuations about the average wind speed over distances of about 
127 meters on the sea surface for large 0. The hypothesis that a large 
part of radar sea return is in practically instantaneous equilibrium with 
the local wind is substantiaed by this analysis. The result of the analy- 
sis of these three time (or space) series are perfectly predictable from 
Tables II and III. Wind variability is large enough to be at least as impor- 
tant as fading at 35°, whereas at 15°, fading completely submerges any 
wind variability. 

These same time series techniques can also be applied to the series 
of original values of or° at 10°, 15", 25° and 35°. The spectra and co- 
variance functions are remarkably like the ones for the transformed wind 
speeds despite the exponential transformation. From these spectra by 
means of equation 7, the effective number of degrees of freedom in the 
sample of 189 values of or° can be found. 
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It now becomes possible to put one standard deviation confidence 
intervals on the estimates of the mean values of such quantities as 


10 1o §10 


7 ° 05 °) 

cr ° (10 •) 


as given in Table IV. The effective number of independent points instead 
of being 189 has been reduced to the values shown, but the loss of de- 
grees of freedom is much less than in the anemometer record in which 
2000 data points were equivalent to 178 independent points . It has been 
assumed that cr o (l0 o ) is independent of the other cr 0 values so the vari- 
ances have been added and the lowest applicable degrees of freedom has 
been used. 


TABLE IV. - ANALYSIS OF FLIGHT 7 P RUN 1, LINE 1. 


1 89 data points 


0 



10 

1 — ° 

log „ O' 
s 10 

SD 

( 

r (db) 

VAR 

N 

\/N* 

10 

o 



10.28 

0. 

823 

0.68 

152 

12 

15 

o 



9. 100 

0. 

807 

0.65 

143 

12 

25 

o 



1.596 

0. 

975 

0.95 

142 

12 

35 

o 



-7.056 

1 . 

081 

1.18 

95 

10 

35 

o 

10 

o 

-17. 34 

1 . 

36 

1 . 86 



25 

e 

10 

o 

- 8.68 

1 . 

28 

1 . 63 



15 

o 

10 

0 

-1.18 

1 . 

15 

1 . 33 



10 

iogj 

1 

0 

( 7(35 *) / 

7 (io“) 

) = 

-17. 34 ± 

0. 136 


10 

logj 

l 

0 

( 7(25°) / 

7(io°) 

) = 

-8.68 ± 0 

. 107 


10 

iogj 

1 

0 

( 7(15“) / 

7(io“) 

) - 

k— * 

00 

1+ 

o 

.096 



The normalized average c values for 35 °, 25 °, and 15 ° have been 
e stimated to within +.0,136, _+ 0 . 1 0 7, and _+ 0 . 096 db r e spectively . 

It will be recalled that the mean of the anemometer record was 
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17 knots the standard deviation was 1.95 knots. The degrees of freedom 
were estimated by one procedure to be 178 for an 8000 second (133 
minute) record. The flight line was shown to be the equivalent of a 49 
minute anemometer record so that the sampling variability in the esti- 
mate of the mean would have (49/133) X 178, or about 65 degrees of 
freedom. This results in an estimate_of the one standard deviation con- 
fidence interval on the mean wind of v = 1 7 ± 0. 24 knots. 

This estimate of the variability of the mean wind is suspect because 
of the shape of the wind spectrum. It was checked by a further study of 
the anemometer record at Argus Island. Actual 49 minute overlapping 
samples of three portions of the record were averaged for times near 
Flights 9 and 7. The actual standard deviations of these averages were 
0. 33 knots, 0. 38 knots and 0. 55 knots. The range in the averages was 
1 . 24 knots, 1 . 41 knots and 2. 25 knots. The values of 0. 38 knots and 
1.41 knots apply to Flight 7, Line 1, Run 1. 

Similar results have been given by Bradley (1971) for other data 
runs as well as this one. The method he used to estimate the standard 
deviation of the sample mean of c r° was more conservative so that 
broader confidence intervals were obtained. The estimated variability 
of the winds is comparable. Nevertheless, even these confidence inter- 
vals are quite narrow. 

All of these results have to be corrected for atmospheric stability 
and the elevation of the anemometer and for variations in the antenna 
pattern both within missions and from one mission to another. 

For 9 = 35 °, the results are summarized for most of the data ob- 
tained on Missions 119 and 156, after all appropriate corrections, in 
Table V. There are 20 data sets identified by flight and mission. The 
meteorologist's estimate of the wind is given by V]y[, good to _+ 0.5 knots. 
The value of cr^ is also given along with the wind computed from this 
radar measurement, V R> as tabulated. The difference to one -tenth 
of a knot is also tabulated, and the root mean square value of the differ- 
ence in 2.6 knots. 

A scatter plot of versus is shown in Figure 14. The scatter 

is really not very large and is due to both the difficulties of specifying 
the wind from the meteorological point of view and the scattering cross 
section from a radar point of view. 
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V = Meteorologically Determined Wind 

M 



BY RADAR. 
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TABLE V.- METEOROLOGICAL WIND ESTIMATES, cr° VALUES AND 
WIND SPEED COMPUTED FROM cr° VALUES FOR UP- 
WIND 35 0 INCIDENCE ANGLE FOR DATA FROM MISSIONS 
119 AND 156. 


Set 

Mis sion 

Flight 

V 

o- °(db) 

V R 

Differences 

1 

119 

9-3-3 

6 

-23. 1 7 

5. 7 

0. 3 

2 

119 

9-2-2* 

6 

-21 . 23 

7.6 

-1.6 

3 

156 

3-4-5 

10 

-20.02 

9.2 

0.8 

4 

156 

3-4-6* 

10 

-21 . 04 

7.9 

2. 1 

5 

119 

5-1-2 

12 

-18.02 

12. 6 

-0.6 

6 

156 

5-4-5 

15 

-16.65 

15.5 

-0.5 

7 

156 

5-4-6* 

15 

-17. 47 

13. 7 

1 . 3 

8 

119 

7-1-1 

16 

-16.25 

16.5 

-0.5 

9 

119 

7-1-2* 

16 

-15.01 

20.0 

4 

10 

119 

9-1-1 

16 

-16.88 

15.0 

1 

11 

119 

9-1-2* 

16 

-15.5 

18.5 

-2. 5 

12 

119 

7-1-7 

17 

-15.70 

17. 9 

-0.9 

13 

119 

9-1-3 

17 

-16.85 

15.0 

2 

14 

156 

7-8-5 

22 

-13.67 

24.5 

-2.5 

15 

156 

■^1 

1 

00 

O ' 

22 

-14.97 

20.0 

2 

16 

119 

2-1-1 

22 

-15.69 

18.0 

4 

17 

119 

2-1-2* 

22 

-15.52 

18.5 

3.5 

18 

156 

2-8-1 

31 

-11.23 

35. 7 

-4. 7 

19 

119 

3-1-2* 

32 

-10.96 

37. 3 

-5. 3 

20 

119 

3-1-1 

33 

-12. 15 

31 . 0 

2 


RMS ERROR 2.6 

knots 


An Inc onsistency . - The data on Flight 7-1-1 show that the scatter 
in the measurement of a particular average scattering cross section is 
small. The scatter in the measurement of the average wind is also 
small. Moreover, the data show that the measurements of c r° are 
actually responding to the gusts and lulls in the wind over the ocean sur- 
face along the path scanned by the beam. 

On the other hand, the data in Table V and Figure 14 show a much 
larger scatter between the sea return and the wind, or the wind com- 
puted from the sea return and the meteorological wind, than can be ex- 
plained on the basis of measurements made on a particular flight line 
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and run. Even if these estimates are off so that <r 0 scatters by _+ 0. 3 db 
and 17 by ± 1 knot, the scatter in the full data set is still too large. 

The greater scatter in the total set of data with rms errors of 0. 7 db 
and 2.6 knots, or so, must be explainable in terms of factors other than 
the duration of the flight line and the concurrent measurement of the wind, 
especially for the Argus Island, Mission 119 results. These sources of 
scatter can be instrumental and meteorological and, or, oceanographic. 

For example, the starred values for the flights in Table V were ob- 
tained with the aircraft flying downwind looking upwind with the aft point- 
ing beam. Many of these values go with a corresponding set of values 
with the aircraft flying upwind looking upwind with the forward pointing 
beam. These pairs of measurements were made within just a few min- 
utes of each other. Although Bradley (1971) attempted to calibrate out 
fore beam-aft beam antenna pattern differences, these values suggest 
the possibility of a residual error in this effect. 

In processing the scatterometry data, one important part of the cal- 
culation is the ground speed of the aircraft, which is used in determin- 
ing the Doppler shift. This in turn is used to calculate the backscatter 
for each angle. A study of the determination of the ground speed for 
Mission 1 19 lead to the conclusion that the ground speed was not deter- 
mined with complete accuracy. Revised values of the ground speed are 
given by Pierson and Cardone (1971). The use of an incorrect value 
for the ground speed has two effects. It shifts the measured value to 
a value of 9 different from the one desired, and it gives a value that 
needs to be corrected for by a different angle in the antenna gain cor- 
rection. The combined effect is rather difficult to account for fully, 
but it can surely introduce additional scatter. 

Other sources of residual error can be meteorological and/or 
oceanographic. The form of the wind profile with height and the drag 
coefficient must be known to be able to compute the variation of wind 
with height and the friction velocity. The scatter in such a basic quan- 
tity as the drag coefficient has been documented for Argus Island ob- 
servations by DeL-eonibus (1971). 

There is also the problem of the effect of the Argus Island Platform 
and supporting structure on the wind flow past the anemometer on the 
tower. This effect was studied by Thornthwaite et al (1965) at Argus 
Island and by Mollo-Christenson and Seesholtz (1967) for a similar. 
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but smaller tower in a wind tunnel. For a south wind, the speed mea- 
sured by the anemometer could be 4 or 5$ too high. However, not enough 
information is available to make corrections for tower effects for the 
variety of wind directions encountered during Mission 119. Hence, no 
correction for this effect was attempted. The variation of this effect 
with wind direction could easily account for another ^ to 1 knot scatter 
in the estimate of the average wind. 

There is one oceanographic source of scatter that has not been re- 
moved from the data. The high frequency part of the wave spectrum is 
in equilibrium with the local wind but the lower frequency gravity wave 
portion need not be in equilibrium, as a fully developed sea, with the 
wind. Preliminary models of radar sea return show a dependence of <r° 
on the longer higher waves in the spectrum and the present relationship 
can eventually be corrected for this variation in more complicated systems 
by means of both theory and multiple regression techniques. At the 
present level of development of the theory of sea return this added refine- 
ment is probably premature. 

It is our opinion that much of this scatter has to do with meteorologi- 
cal, aircraft and radar instrument difficulties. With scanning pencil 
beam radars such as the one built under the AAFE Radscat program and 
for SI 93 on Skylab, a good part of this scatter may not be present. Under 
the AAFE Langley Radscat program, plans are being made to minimize 
the problems of estimating the average wind. 

To study this problem further, the two sources of variability identified 
in the preceding sections will be extrapolated to spacecraft altitudes for 
SI 93. The result will be predictions of the capabilities of SI 93 in clear 
air to determine the wind speed and direction over the oceans under the 
assumption that the scatter in the data will be caused by the combination 
of fading and a turbulent wind over the water. 


EXTRAPOLATION TO SI 93 


The two sources of scatter in the measurement of <r° obtained with 
aircraft that were identified in this report have their counterpart for 
SI 93 and the effects of both fading and turbulence can be taken into ac- 
count for measurements to be made by Si 93. Both of these effects re- 
sult in a very small scatter in the measurements. 
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Skylab SI 93 Scanning Patterns . - Three typical orbit segments for 
Skylab during which Si 93 might be turned on are illustrated in Figure 15. 
The wind field near the surface of the ocean as deduced from ship reports 
is illustrated in Figure 16 and the cloud cover that was present is shown 
in Figure 1 7. 

Orbit segment A in Figure 15 shows the cross-track non-contiguous 
scanning mode. The antenna points straight down and at angles of 15.6°, 
29. 4°, 40. 1 and 48 0 off of the vertical at the spacecraft. Because of the 
curvature of the earth the appropriate 0 values will be larger than these 
angles . 

The illuminated areas shown by the open circles are not drawn to 

scale. They are larger as the angle increases. At each cell, six differ - 

ooo o ^ , 

ent measurements will be made; Oyy, cr^^, anc * ^VH anc * ™ le P asslve 

microwave temperature for two polarizations. The pass shown would 
obtain 8 70 separate measurements in about seven minutes. 

Orbit segment B shows a northbound pass and orbit segment C 
shows a southbound pass on the non-contiguous forward scanning mode. 
Each black ellipse represents a series of six measurements at each of 
the five different angles at nearly the same place on the ocean surface 
for a total of 30 measurements. 

Figure 16 shows a wind field analysis such as might occur during a 
Skylab pass. The various scanning modes would scan cells on the ocean 
surface in which the wind speed varied all the way from 5 knots, or so, 
to perhaps 50 or 70 knots, depending on the scanning mode and the exact 
location of the radar illuminated areas on the sea surface. 

The M radar scatter 11 or ’'fading variability 11 of SI 93 . - S193 measure- 
ments have a source of scatter analogous to the fading that results from 
the 0. 3275 second averages in the aircraft data. Each of the four SI 93 
radar measurements are averaged over about 0.6 seconds and the band- 
width of the system is such that the received power is measured to within 
7$ for each measurement. Thus within one standard deviation the radar 
scattering cross section is good to within 1 + 0.07 of the expected value. 

In fact, the averaging time is a function of the antenna angle and has been 
chosen so that this value is an upper bound for all scan angles. 
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FIG. 15. CROSS TRACK AND ALONG TRACK SCANNING PATTERNS 
FOR S193 ON SKYLAB. 




6 . 
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o This design fact yields the standard deviation of measurements of 
o- which is ± 0. 3 db. Table VI shows the results using this value for 
SDg. in equation 13 for the off- vertical values of 0 and the values of B 
found by Bradley (1971). The upwind B for 15° was used for 15. 6 
the B for 25° was used for 29.4° and the B for 35 “ was used for both 
40. 1 * and 48°. In actuality, the values for B for these angles are likely 
to be higher, which in turn implies a lower percentage variability in the 
wind estimates. 

It also must be assumed that the same kind of equation as equation 1 1 
will also be valid for non-normalized measurements of cr° at these vari- 
ous values of 9. Since the transmitted power is monitored in Si 93, the 
values of cr° should be primary measurements that will not have to be 
normalized to 10 ® as in the aircraft flights. The constant. A, will be 
different for the SI 93 measurements but the equation should still be of the 
same form. Results similar to Table VI can also be computed for cross 
wind conditions and somewhat larger percentages of scatter would result. 

TABLE VI.- CONSERVATIVE ESTIMATES OF PER POINT 
VARIABILITY DUE TO FADING OF UPWIND 
WIND SPEED MEASUREMENTS BASED ON 
SINGLE MEASUREMENT OF tr° BY SI 93. 


log 1.07 = 0.030 

SD =0.3 db 
tr 


0 

l°g l0 (1 + Av/v) 

1 + A V/V 

15.6° 

0. 0909 

1 . 23 

29. 4° 

0.0266 

1 . 06 

40.1* 

0.0201 

1 . 05 

o 

00 

0.0201 

1.05 


In essence for angles of 29.4° and greater, the wind speed estimate 
for a single measurement of oyy will scatter by about ± 6$ about the true 
value due to this source of variability. All theoretical and observational 
data to date suggest that ^j-jv an< ^ °'VH W H1 a-11 provide additional 

equally independent estimates of the wind-speed and that the scatter will 
be comparable. The average of four such estimates, since each one will 
be independent, will produce a final estimate within 3$ of the expected value. 
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The available data on <Tyy (0) as a function of wind speed all suggest 
that at 9 = 0, jyy(O) decreases with increasing wind speed and that at an 
angle near 10 % or so, 0yy (0) does not depend strongly on wind speed. 

This is why the normalized measurements at 15° are relatively insensi- 
tive to wind speed and why the values for 1 + AV/V and 1 - AV/V are 
large for 15 An operational instrument could avoid this problem by 
tilting the antenna forward as described by Moore and Pierson (1971). 

The values given for 15* in Tables I, II, III, IV, and VI are all a natural 
consequence of the shape of the scattering cross section curve and are 
therefore to be expected. 

S193 also measures passive microwave temperature in two polariza- 
tions. Moore et al (1971) have combined the data obtained by Nordberg 
and Hollinger in a plot that shows that, given the sea surface temperature, 
AT b ( °K), increases at 0. 96 + 0. 16 °K/(m/s) for 19 GHzas the wind varies 
from zero to 30 meters per second. 

The variation of AT bhh (°K) for 1 3. 9 GHz will not be as strong as 
that at 19 GHz, but througn cloudless skies after correction for atmos- 
pheric temperature and water vapor and without obvious corrections for 
fetch and duration, which change the amount of foam, a similar dependence 
about 75$ as strong could lead to the dependence given by equation 16 

(0.78i 0.16)AT BHH (*K) (16) 


It may then be possible to estimate the wind speed at 1 5 * to within about 
20$ by a passive measurement. Refinements in calculating the amount 
of foam could reduce this scatter substantially. At 0°, there will again 
be a wind speed dependence for cr ® which can be combined with the pas- 
sive microwave to determine the wind. 

Moreover, the passive microwave measurements appear to be inde- 
pendent of wind direction. It should be possible to devise a proceedure 
in which all six measurements for a 0 greater than 15 ° can be used to 
determine both wind speed and direction to within a choice of four possible 
angles, + X and +X + 180°, at each cell through cloudless skies. Three 
of these four angles can often be eliminated by synoptic considerations. 
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— -fepts °f Turbulent Wind Field on Measurements Made with Si 93 . - 
The patch of ocean illuminated by the radar beam on Si 93 varies in area 
and dimensions as a function of the angle of the beam off the vertical. 

The area of the patch is 95 km 2 for 0° and increases to 385 km 2 at 48 * 
as shown in Table VII. 

The area given in this table is computed on the as sumption that the 
main lobe of the antenna pattern illuminates the sea surface with constant 
power over the entire given area and that there is zero power outside 
the area. It is also an instantaneous value. The time required to make 
the measurement spreads out the illuminated area and applies a trape- 
zoidal weight function to a broadened pattern. This incremental area is 
also tabulated for the cross track mode, but it was not used in the rest 
of the calculations. This incremental area actually increases the value 
of N by an amount equivalent to adding about half the incremental area 
to the values tabulated under "area". The rest of the entries in the table 
are therefore conservative. 

The analysis of the Argus Island anemometer record led to the con- 
clusion that for a wind speed near ]6 knots the anemometer record would 
be the equivalent of measurements of wind speed every 350 meters along 
a line on the ocean surface in the wind direction. Lumley and Panofsky 
(1964) show that the horizontal wind decorrelates in the cross wind direc- 
tion at least as fast as in the upwind -downwind direction so that points on 
a square grid 350 meters on a side over the area illuminated by the radar 
beam could represent the equivalent number of independent observations 
of the wind. The wind at each of these points, and in fact over the entire 
area, is fluctuating above and below the desired average value as repre- 
sented by equation 5 in much the same way that anemometer records 
fluctuate. 

A square 350 meters on a side has an area of (0. 35) 2 = 0.123 square 
kilometers so that there are effectively 8 independent points per square 
kilometer of illuminated area for a 16 knot wind. 

The number of effectively independent points for a 16 knot wind is 
shown in Table VII for each antenna angle. These values have to be re- 
duced by a factor of four when the wind speed is doubled and by a factor 
of nine when the wind is tripled because of the scaling laws for turbulence. 

Even for a 48 knot wind, there are at least 84 degrees of freedom in the 
vertical. 
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TABLE VII. - CONSERVATIVE DEGREES OF FREEDOM (N) IN 
ESTIMATES OF VARIANCE OF THE AVERAGE 
WIND OVER THE RADAR ILLUMINATED^ PATCH, 
AND £5$ CONFIDENCE INTERVAL ON V FOR 


SD /V 

V 

= 0. 20 AS GIVEN BY 

V + (2SD 

V 

A/N). 


Angle 

0 ° 

15.6* 

29. 4° 

40.1° 

e 

00 

Area (KM) 

95 

108 

151 

236 

385 

Incremental 






Area (KM)^ 

21 

39 

68 

97 

142 

N (V = 16 knots) 

760 

864 

1208 

1888 

3080 

N (V = 32 knots) 

190 

216 

302 

472 

770 

N (V = 48 knots) 

84 

96 

134 

209 

347 

1 / \/N (16) 

. 036 

. 033 

. 028 

. 023 

. 018 

1 /-s/N (32) 

.072 

. 066 

. 056 

. 046 

.036 

1 /\fN (48) 

. 109 

. 102 

. 086 

. 069 

. 054 


95$ Confidence Interval 


1 6 knots 

16.24 

16.20 

16.18 

16. 14 

. 16.11 


15. 76 

15.80 

15.82 

15.86 

15.89 

32 knots 

32. 92 

32.85 

32. 72 

32.59 

32. 46 


31 . 08 

31 . 15 

31 . 28 

31 . 41 

31 . 54 

48 knots 

50.09 

49. 96 

49.65 

49. 32 

49. 03 


45.91 

46.04 

46.35 

46 . 68 

46. 97 


If the standard deviation of the wind is 20$ of the mean wind, 95$ 
confidence intervals on the mean wind can be found by computing as in 
Table VII. 

These confidence intervals are delightfully narrow. For inter- 
comparison with other numbers in this paper, they should be halved. 
The 95$ confidence intervals are plus or minus less than 1 knot for all 
entries in the table except 48 knots. The largest percentage variation 
for 0 0 and a 48 knot wind is 4. 4$. It can be concluded that the illumi- 
nated patch is large enough to provide a truly stable estimate of the 
average wind for the purposes of synoptic meteorology. 
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The dimensions of the spot illuminated by the radar beam are given 
in Table VIII. The spot at 48° is roughly elliptical with the doubled 
major axis distance equal to 29 km and the doubled minor axis distance 
equal to ] 7 km as an example from the table. 

Here again the spot is actually stretched in one direction due to 
spacecraft motion. The stretching makes the values in the table con- 
servative . 

TABLE VIII. DIMENSIONS OF PATCH OF OCEAN ILLUMINATED 
BY SI 93 RADAk BEAM AND FREQUENCIES IN AN 
ANEMOMETER RECORD THAT CORRESPOND TO 
THESE DIMENSIONS. 


Angle 
Area (KM) 

Approximate Length (KM) 
Approximate Width (KM) 

0 ° 
95 
1 1 
1 1 

15.6 ° 
108 
12 
12 

29. 4° 
151 
15 
1 3 

40. 1 ° 
236 
21 
15 

48 ° 
385 
29 
17 

Equivalent Time Minutes 
8 meters/sec 

23 

25 

31 

44 

60 


23 

25 

27 

31 

35 

16 meters/sec 

11 

12 

15 

22 

30 


11 

12 

1 3 

15 

17 

24 meter s /sec 

7 

8 

10 

15 

20 


7 

8 

9 

10 

12 

Approximate Cycles /Hour 
8 meters /sec 

2 . 6 

2. 4 

1 . 9 

1 . 4 

1 . 0 


2.6 

2.4 

2. 2 

1 . 9 

1 . 7 

1 6 meter s / sec 

5.5 

5.0 

4. 0 

2. 7 

2. 0 


5.5 

5.0 

4.6 

4.0 

3.5 

24 meter s /sec 

8.6 

7. 5 

6.0 

4. 0 

3. 0 


8.6 

7. 5 

6. 7 

6.0 

5.0 
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Consider a sinusoidal perturbation in the wind with a length in space 
given by these tabulated dimensions and assume that it is advected by the 
mean wind. Such a perturbation moving at 8 meters per second would 
take 60 minutes, as tabulated, to travel 29 km across the beam. The 
29 km ocean distance is thus the equivalent of a perturbation with a 
period of 60 minutes in an anemomber record when the mean wind is 
8 meters/sec. 

The periods that result for each antenna angle and for the two dis- 
tances across each patch are all given in Table VIII. They range from 
7 minutes to 60 minutes. The reciprocals of these numbers converted 
to cycles per hour are given as the last set of entries. The values range 
from 1.0 to 8.6 cycles per hour. Either way, when these periods, or 
frequencies, are located on the spectrum in Figure 3, they all fall in 
the valley between the microscale turbulent fluctuations and the synoptic 
scale winds. The illuminated patch of sea surface from the scanning 
pencil beam on SI 93 acts as a very efficient way to obtain the quantity 
shown in equation 5 by acting as a low pass filter to suppress the micro- 
scale fluctuations and provide a very stable scattering cross section 
measurement representative of the desired synoptic scale wind 1 . 

The two dimensional "box-car", or square sided, weight function 
assumed for the antenna pattern is an over simplification of the problem. 
The power illuminating the surface decreases in a radial direction from 
some center spot whose location depends on geometry and the movement 
of the spacecraft during the time required for one measurement. To ac- 
count for 96 *j> of the total power involved, the dimensions of instantaneous 
pattern must be increased by a factor of 1.65 so that the area is increased 
by a factor of about 2. 6. The exact form of the weight function will de- 
pend on many factors. It will, however, be a smoothly varying function 
that builds from zero to a unimodel peak, perhaps with roughly elliptical 
lines of equal power. 

Such a weight function when used in equation 5 would produce a filter 
in vector wave number space that would act like a fairly sharp low pass 
filter and provide a good measurement of the synoptic scale wind. Work 
on the antenna pattern is in progress at the University of Kansas, and 
it will be possible to refine these calculations for application during 
Skylab . 
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In summary, two of the sources of variability in estimates of the 
mean wind for synoptic scale purposes, namely fading in the scattered 
signal and turbulent fluctuations in the wind over the water are both ade- 
quately averaged in the design of SI 93. The two sources of scatter com- 
bined produce a total scatter of, at most, 7$, except at 15°, so that two- 
thirds of the time a 30 knot wind would be estimated to be between 28 and 
32 knots and a 40 knot wind would be estimated to be between 37 and 43 
knots . 

Additional Sources of Scatter for S193 Data. - The root mean square 
error for <r° in db found by Bradley (3 971 ) as given in Table I, and the 
scatter in Table V and Figure ]4 was not all explained by fading effects 
and the turbulence in the wind. For both SI 93 and the airborne AAFE 
Langley Radscat, the instrumental problems due to the use of the fan 
beam Doppler system will not be present. There remains the problem 
of determining the wind velocity accurately for the areas of the ocean to 
be sampled by these instruments. 

Unfortunately, the improved design of these two new instruments 
cannot be matched by an improvement in the accuracy with which the 
wind speed can be determined meteorologically, especially for SI 93. 

Some improvement can be obtained for experiments with the airborne 
instrument. 

The problem of the adequate determination of the synoptic scale wind 
for patterns scanned by Si 93 as in Figure 15 has been treated in a pro- 
posal to NASA by McClain, Moore, Pierson and Talwani (1971) and by 
an amended proposal in which a concurrent aircraft program is described. 
The wind field as shown in Figure 16 ought to be smoothly varying over 
synoptic scale distances. Computer based procedures such as those 
developed by Cardone (1969) and described by Pierson (1970) can be used 
to smooth the wind field given by ship reports and surface atmospheric 
measurements so as to determine the wind. It is to be expected however 
that these estimates of the wind velocity will introduce considerable 
scatter in the regressions of cr° values on wind velocity. Not only will 
the upwind and crosswind situations have to be treated, but also all pos- 
sible angles between the beam and the wind direction. The data base for 
these studies will, however, be much greater so that the scatter will be 
more easily understood. 

The use of SI 93 on Skylab has converted an instrument system, which 
when used on an aircraft, would be responsive to microscale effects and 
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therefore have concurrent measurement and sampling problems, to a 
true synoptic scale instrument in which the microscale is averaged out. 
The result may even be that winds determined by Si 93 will be so good 
that other wind measurement systems such as ship reports and anemo- 
meters on buoys will be calibrated against SI 93 instead of the other way 
around. 

Other Problems of Si 93. - The preceding considerations have all 
been carried out under the assumption that the beam passes through clear 
air to the sea surface. There is the additional problem of determining 
the wind velocity through clouds and precipitation as illustrated by Figure 
17. This problem has to be treated by Moore et al (1971). The essential 
points are that passive measurements lose contact with the sea surface 
when clouds are present long before the radar measurements and that the 
passive measurements can be used to correct the radar measurements for 
atmospheric cloud effects and find the wind at the surface. 

Also there is the possibility that variations in radar backscatter can 
be caused by variations in sea surface temperature, salinity, and surface 
properties which produce changes in the dielectric properties of the water. 
These effects may have to be calibrated out. Finally, theories of radar 
sea return show that it depends on the entire wave spectrum and not on 
just the high wave number part of the spectrum. The dominant effect is 
from the short waves, but corrections for that part of the spectrum not 
in local equilibrium with the wind will have to be developed. These cor- 
rections are expected to be small, but necessary. 


CONCLUSIONS 


The extrapolation of laboratory and aircraft radar sea return data to 
spacecraft altitudes has shown, for Si 93 on Skylab, that fading effects 
will contribute very little scatter to the measurements and that the micro- 
scale turbulent fluctuations in the wind will be averaged out. The result 
should be very stable estimates of the synoptic scale winds. 

Additional sources of scatter in the aircraft program were identified 
as being instrumental, associated with the fan beam Doppler radar that 
was used and the anemometer location on Argus Island, and meteorologi- 
cal and oceanographic. The radar problems will not be present for SI 93. 
There will still be problems in defining the wind by means of meteorologi- 
cal measurements to comparable accuracies. 
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Other problems, unique to a space instrument, such as making mea- 
surements through clouds and precipitation were identified. Methods for 
their solution were given. 

The volume, and the improved quality, of the data from SI 93 should 
make it possible to meet the scientific objective for Si 93 which is to prove 
conclusively that simultaneous measured values of the radar scattering 
cross section matrix and the microwave temperature for two different 
polarizations at 13.9 GHz will provide data from which the winds in the 
planetary boundary layer over the oceans can be inferred. 
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